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We use Floquet theory to study the maximum value of the stroboscopic group velocity in a one¬ 
dimensional tight-binding model subjected to an on-site staggered potential varying sinusoidally in 
time. The results obtained by numerically diagonalizing the Floquet operator are analyzed using 
a variety of analytical schemes. In the low frequency limit we use adiabatic theory, while in the 
high frequency limit the Magnus expansion of the Floquet Hamiltonian turns out to be appropriate. 
When the magnitude of the staggered potential is much greater or much less than the hopping, 
we use degenerate Floquet perturbation theory; we find that dynamical localization occurs in the 
former case when the maximum group velocity vanishes. Finally, starting from an “engineered” 
initial state where the particles (taken to be hard core bosons) are localized in one part of the chain, 
we demonstrate that the existence of a maximum stroboscopic group velocity manifests in a light 
cone like spreading of the particles in real space. 
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I. INTRODUCTION 

In the real world, the speed at which information can 
propagate is limited by the speed of light; this results 
in the light cone effect as postulated by the special the¬ 
ory of relativity. Is there a similar upper bound of the 
speed at which correlations (information) can propagate 
in interacting quantum many-body systems? Following 
the seminal work by Lieb and Robinson^, which estab¬ 
lished the existence of a maximum group velocity in a 
one-dimensional spin chain with a finite range interac¬ 
tion, some recent studies have explored this conjecture 
in several interacting many-body systems; these studies 
do indeed exhibit an effective light cone that sets a bound 
on the speed of propagation of correlations. This is re¬ 
flected for example, in the growth of block entanglement 
entropy following a quench^, or the collapse and revival 
of the Loschmidt echo^. The light cone like propagation 
of quantum correlations has also been observed experi¬ 
mentally by quenching a one-dimensional quantum gas 
in an optical lattice^. 

In parallel, there have been a plethora of studies of 
closed quantum systems driven periodically in time in 
the context of defect productions^ii, dynamical freezing^, 
dynamical saturation^ and localizationiSrJ^, dynamical 
fidelityi^, and thermalizationi^ (for a review see Ref. [H). 
The study of periodically perturbed many-body systems 
has also gained importance because of the proposal of 
Floquet (irradiated) graphenei^^— , Floquet topological 
insulators and the generation of topologically protected 
edge statesiSr.^ some of which have been experimentally 
studied^^— . 

In this work, we use Floquet theory to explore the stro¬ 
boscopic (i.e., measured at the end of each complete pe¬ 
riod) group velocity of a system of hard core bosons resid¬ 
ing on a one-dimensional lattice in the presence of a stag¬ 
gered potential which is varying sinusoidally in tim o^^d^ . 


In particular we study the maximum value of the group 
velocity to observe the consequent light cone effect. Al¬ 
though the time-independent version of the model is in- 
tegrable, the periodic sinusoidal perturbation renders the 
situation rather complicated since the corresponding Flo¬ 
quet operator cannot be obtained in a closed analytical 
form unlike the case of periodic perturbations which are 
piece-wise continuous in tim o^^'^^ . One therefore has to 
use various approximation schemes valid in the appropri¬ 
ate regions of the parameter space to analyze the behav¬ 
ior of the stroboscopic group velocity. 

The paper is organized in the following way. In Sec. m 
we present the Hamiltonian of the model under consider¬ 
ation and discuss the generic behavior of the maximum 
value of the stroboscopic group velocity 14nax as a func¬ 
tion of the amplitude and frequency of the periodic per¬ 
turbation; this is derived using the numerically obtained 
Floquet quasi-energies. In Sec. imi we use adiabatic the¬ 
ory to find the behavior of Vjnax in the low frequency limit 
while the high frequency limit is treated within a Mag¬ 
nus expansion in Sec. [TV] In Sec. |V](Sec. IVII) . we use a 
Floquet perturbation theory^ when the hopping term is 
much smaller (greater) than the amplitude of the stag¬ 
gered potential. In the case of a large magnitude of the 
staggered potential, we point to the situations when the 
maximum stroboscopic group velocity vanishes resulting 
in the so-called dynamical localization. We demonstrate 
a light cone like propagation of particles in real space in 
Sec. IVIII Concluding remarks are presented in Sec. IVIIII 

II. MODEL AND THE STROBOSCOPIC 
GROUP VELOCITY 

We consider a Hamiltonian of the tight-binding form 
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FIG. 1 : (Color online) Density plots showing the variation of the maximum value of the stroboscopic group velocity as a 
function of the driving frequency cj and the magnitude of the sinusoidally varying staggered potential Vo (both in units of 
7 = 1 ). In (a) V)nax is always finite (red region) for large u and small Vq (top left corner). But there are some regions where 
Vmax becomes very small (blue region) in the intermediate cj range and for large Vo (bottom right corner), (b) shows that there 
exists an array of regions where I 4 nax vanishes (dark blue regions). The line-like regions of zero velocity are not equispaced 
in to for a given value of Vo. (c) shows that Fmax never becomes zero for small u) and small Vo, although Fmax is very small 
for extremely low frequencies (blue region near the horizontal axis), (d) shows that t^nax decreases as one increases Vo while 
keeping uj fixed at small values. At extremely small frequencies, Knax shows behavior similar to that in (c). (e) A schematic 
diagram showing the validity of the various limits of the Floquet theory in different regions in the uo — Vo plane. 


where 7 is the hopping amplitude, and 6; denotes bosonic 
annihilation operators defined on a one-dimensional lat¬ 
tice satisfying the hard core condition (fej)^ = (6;)^ = 0. 
The single particle dispersion is given by Ek = —27 cos k. 
This model describes a system of hard core bosons in 
a gapless superfiuid phased. (By the Jordan-Wigner 
transformation^, this system is equivalent to one with 
spinless fermions). When the system is perturbed by a 
spatially alternating potential V) = (—l)*Vo, a gap opens 
up in the spectrum for any non-zero value of Vo thereby 
driving it to a gapped Mott insulator phase. Our aim 
here is to investigate the response of the system sub¬ 
jected to an alternating potential varying sinusoidally in 
time as V{t) = Vosm{ujt). We will analyze the behav¬ 
ior of the stroboscopic group velocity, measured after n 
complete periods of the driving, and its maximum value 
Vjnax as a function of the driving frequency w and the 
amplitude Vq. 

In the presence of the alternating potential, the Hamil¬ 
tonian in Eq. 0 reduces in momentum space to a 2 x 2 
matrix form in terms of the momenta k and k + tt, 

Hk{t) = —2^008 k + Vo sin{u!t) , (2) 

where tr^’^ denote pseudo-spin Pauli matrices. Clearly 
the spectrum is gapless for Vq = 0 . We will set 7 = 1 , 
Planck’s constant h = 1 and the lattice spacing equal to 
1 in the rest of the paper. Hence to, Vo, k and Knax (to 
be defined below) will all be dimensionless. 


Using the Jordan-Wigner transformation, the time- 
independent part of the Hamiltonian can be mapped 
to a system of spinless fermions on a one-dimensional 
lattice with a hopping amplitude 7, while the time- 
dependent part of the Hamiltonian corresponds to a stag¬ 
gered chemical potential which is sinusoidally driven with 
a frequency lo. This equivalence between the hard core 
bosonic and the spinless fermionic models does not hold 
in higher than one dimension. Nevertheless, in higher 
dimensions there are non-interacting fermionic models 
with a sinusoidally driven chemical potential for which 
the behavior of the stroboscopic group velocity and the 
dynamical localization is similar to what is reported here. 

Defining the time period T = 27 r/w, the stro¬ 
boscopic Floquet operator for each momentum mode 
is given by the unitary operator Tk{Vo,T) = 

Texp (^—i dtHk{t)^ , where T denotes the time¬ 
ordering. This operator cannot be computed analyti¬ 
cally for a sinusoidal driving. One can however numer¬ 
ically calculate J-k and find its eigenvalues which take 
the form exp(—I^^T) where /i^ = are the quasi¬ 
energies. The group velocity can be obtained from the 
quasi-energies as 

Vk = d^ik/dk. (3) 

The maximum of Vk as a function of k gives the quantity 
Umax, which is the main object of interest in this paper. 
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The physical interpretation of the stroboscopic group ve¬ 
locity is that if the quantum correlations are measured 
only at the end of each complete period, they would ap¬ 
pear to propagate with a maximum speed Knax- This 
quantity is presented in Fig. [T] as a function of Vq and w. 

Upon inspecting the results presented in Fig. [U one 
finds that Umax tends to saturate at some value for 
large values of lj and small or intermediate values of Uq 
(U o < 1) [see Fig. [T] (a)]. The maximum group velocity 
shows an interesting behavior when both Vq and w be¬ 
come large, as shown in Fig. |T] (b). In this limit, one 
finds that for a given frequency lu, Umax vanishes in reg¬ 
ular intervals of Vq. On the other hand, for a given Vq, 
the zeros of Unax lie in increasing intervals of w. In this 
regime Umax is given by the zeros of a Bessel function as 
we will show below. But when uj is small. Unax never be¬ 
comes zero, although in the limit a; —>■ 0, Unax becomes 
very small irrespective of U; see the bottom regions of 
Figs, d] (c) and (d). The maximum group velocity grad¬ 
ually decreases with Vq if one keeps uj fixed at a lower 
value, as shown in Fig. [T] (d). In subsequent sections, 
we will use different analytical methods to analyze the 
various behaviors described here. 


Finally, we note that we will consider the entire Bril- 
louin zone ranging from k = —7r/2 to k = tt/2. The time- 
dependent part of the Hamiltonian in Eq. m is quantum 
critical and gapless in the thermodynamic limit for the 
modes with k = ±7r/2. The minimum frequency scale of 
the bare tight-binding Hamiltonian (with = 1) is deter¬ 
mined by the system size l/^)j while the maximum 
frequency scale appears for the modes fc = 0, tt and is 
equal to 2 since Ek = —2 cos fc. 


III. ADIABATIC LIMIT OF LOW FREQUENCY 

The behavior of the quasi-energy can be explained in 
the low frequency limit where the adiabatic theory holds. 
We will choose the basis states as the eigenstates of the 
pseudo-spin operator cr^, i.e., (1 0)^ and (0 1)^. In this 
limit, the product of the time period T and the Floquet 
quasi-energies is equal to the dynamical phase ac¬ 
cumulated over a complete time period T; this is given 

by 


j 


f*27r/a. 
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where E{x) is the elliptic integral^. [We note that the 
Berry phase term vanishes in this problem since the 
closed path traced out by the Hamiltonian in ([2]) as t 
goes from zero to T is a line in the x — z pseudo-spin 
space; such a line covers zero solid angle at the origin 
{x,y,z) = (0,0,0).] Interestingly, the behavior of the 
quasi-energy can be qualitatively explained up to certain 
values of Vq. Let us elaborate on this below. 

In the limit U <C 1 and very small w, we can use the 
form of the elliptic functions to reduce Eq. o to 




47r cos fc 




ttVq sec fc 
4w 




( 5 ) 


A quasi-degeneracy occurs when this dynamical phase 
~ rmr. This condition successfully gives the num¬ 
ber of quasi-degenerate points along with values of 
the quasi-degenerate momenta, fc = ± arccos[(mw -I- 
[see Eig. [2](a)]. (We note that the 
critical modes fc = ±7r/2 are always quasi-degenerate.) 

In the intermediate potential range U) ^ 1, the be¬ 
havior of the quasi-energy is again determined by the 
adiabatic evolution of two-level systems. The number of 
quasi-degenerate points is successfully given by = rmr, 


where is given by 

4cos^fc 8cos^ fclog(4Uo) 4Uo 

Vqoj Uow oj 

The group velocity can be obtained from Eqs. m 
and (0); in the former case, we find Vk — ±(2sinfc — 
(U^ secfctanfc)/8 -I- 0{Vq)) while in the latter Vk — 
±(2 sin(2fc)/(U)7r)-|-4sin(2fc) log(4Uo)/(U)7r)). Therefore, 
we conclude that Unax is (nearly) independent of lu for 
small frequencies; see Eig. [3] We find that Umax = 2 for 
very small Uq while for Vq » 1, Unax ~ 1/U) as shown in 
Fig.[4|for small values of lu. In the limit of small Uq and lu 
(w < 1 /L ), the system hardly senses the periodic driving 
and hence Umax is determined by the bare tight-binding 
Hamiltonian; in the limit of large Uq, the Floquet pertur¬ 
bation theory holds (see Sec. El) which explains the I/Uq 
decay of Umax- 



IV. MAGNUS EXPANSION FOR HIGH 
FREQUENCY 

The time evolution operator describing the 
Schrodinger evolution of a quantum system is given 
by Texp[—i J*dt'H(t')] = exp[A(t)]. In the Magnus 
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FIG. 2 : (Color online) Plots showing the variation of the 
quasi-energy as a function of k with the potential for 
Vb < 1 . In this parameter regime we choose two cases: (a) 
Vo > uj and (b) Vo < uj. The behavior observed in (a) can be 
explained using adiabatic theory with small Vo. The locations 
of the quasi-degenerate points in (b) are discussed in Sec. IVII 



the following forn>iS: 

= ~ 4r ^ dt2[H{ti),H{t2)], 

=-^l dt, I dt2l dt, 

{[H{h),[H{t2),H{t3)]] + 1^3). (7) 


As shown below, the nth-order term decreases as 1 /w” ^, 
and is therefore vanishingly small for co ^ oo. For the 
model given by Eq. ([2), = aa^, ~ 2Vba/w 

and = —(da^VoTr/Sw^) — {iaV^/u:'^) . 

We will work in the high frequency limit and retain 
terms up to the order 1/w. We then arrive at an effective 
Hamiltonian given by 


Hfk = (« - 


3a Vk 




2aVn 


da^VoTT 

3w2 


( 8 ) 


where a = — 2cosfc. The effective quasi-energies, ob¬ 
tained by diagonalizing (|S]) and retaining terms of order 
(Vb/w)^, are found to be = ±-\/a^ — 2a^V^joP- ~ 
±a(l — VqIoj'^)-, the quasi-degeneracy points occur at 
k = ± 7 r /2 where a vanishes. In the limit w —^ oo, 
« ±a . Therefore, the maximum group velocity 
Vmax becomes 2 irrespective of Vq. This can also be ex¬ 
plained simply by noting that the periodically varying 
perturbation in the Hamiltonian ([2]) vanishes on average 
in the high frequency limit. Moreover, for smaller values 
of Vb, Vbiiax reaches its saturation value Vbnax = 2 at a 
smaller value of lu as compared to higher values of Vq . 


V. LARGE POTENTIAL COMPARED TO THE 
HOPPING AMPLITUDE: Vo » 7 


FIG. 3 : (Color online) Plot showing that the maximum 
group velocity Vmax is nearly independent of the frequency 
(or slightly decreasing with increasing lj) at low frequency 
for both small and large values of Vb. Vnax decreases as Vo 
increases. The inset shows that Vmax exhibits some wiggles 
with UJ for high values of Vo. 


expansion, the operator A(t) is decomposed in the 
form A(t) = The advantage of using a 

Hamiltonian periodic in time is that the stroboscopic 
unitary operator (i.e., the Floquet operator) can be 
expressed in the form F{T) = exp(—ii7i?T), where Hp 
is the corresponding Floquet Hamiltonian. Thus, for a 
periodically driven system, the Magnus expansion en¬ 
ables us to express the Floquet Hamiltonian in the form 
Hp = \ where the can be expressed in 


Let us now examine the case where the hopping am¬ 
plitude 7 can be treated as a perturbing parameter in 
the Hamiltonian given by Eq. ([2]) • It is useful to consider 
a unitary transformation which shifts the time depen¬ 
dence of the Hamiltonian to the diagonal term, so that 
the transformed Hamiltonian takes the form Hk{t) = 
Vb sin(a;t)cr^ -|- cos ka^. (We now set 7 = 1 as usual). 
The time-dependent Schrodinger equation in the new ba¬ 
sis can then be written as = Hk{t)\4>f{t)). Di¬ 

viding both sides of the equation by Vb and rescaling t 
to tVo, the Schrodinger equation can be rewritten as 

= sin(—t)|<^+(t))-b —(t)), 

#fcW) = - (t)) + ^ W)- (9) 

We will set w/Vq = a and 2cosfc/Vo = 6 in subsequent 
calculations. The solutions in the zeroth order of w are 
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given by 


= 

4(0) e 


( Cfc (0) e 


i cos(at )/a 
0 
0 

—i cos{at )/ 



( 10 ) 


where c^( 0 ) denote the probability amplitudes of the 
states 10^) at time t = 0. We find that (T)) = (0)) 

and l^^(T)) = |^^( 0 )), implying that these solutions are 
degenerate in Floquet theory. We therefore employ a 
degenerate perturbation theory to include the hopping 
term perturbatively and find the time-dependent coeffi¬ 
cients c^(t), which satisfy the evolution equations 


*C+(t) = bclit) e-*2cos(at)/a^ 

icl(t) = bctit) ( 11 ) 


To incorporate the correction up to first order in the hop¬ 
ping, we substitute c^{t) appearing on the right sides of 
Eqs. (fTTl) by c^(0), respectively. The solution at time 
t = T = 27r/a is then given by 

c+(T) = c+(0) - (*2^6c^(0)/a) Jo(2FoM, 

= c^(0)-(*27r6c+(0)/a) Jo(2EoM. (12) 


Up to first order in the hopping 7, the Floquet operator 
is given by 


. . _ r 1 -{i2TTb/a)Jo{2Vo/uj) 

[_(i27r6/a)Jo(2Vb/a;) 1 J’ 

(13) 

so that (c+(r) c-{T)r = Fk{T)ic+{0) c^{0)f = 
exp(z0^)(c^(O) c,;, (0))^, where (.)^ denotes trans¬ 

pose. Diagonalizing matrix (USD, one obtains 
the eigenvalues , and hence 6^ = = 

±(27r5/a) Jo(2Vb/w) = 2 cosfcJo(2Uo/w)- 

For large Vq and larg e w, the quasi-energy is given by 
~ ±2 cos kyJu/irVo cos(2Vo/u; — 7r/4). The maximum 
group velocity is t4iax = 2^/u^pKVQCOs{2Vo/uJ — 7r/4) 
which that vanishes at 2Vo/a; = (n-|-3/4)7r. This matches 
the observed numerical results presented in Fig. 21 We 
note that the maximum group velocity vanishes when 
Jo(2Uo/u;) = 0. 

Furthermore, given = ±2cosfcJo(2Vb/w), one 
can find the values of the momenta for which quasi¬ 
degeneracies occur in the Floquet spectrum given by 
/ifcT = mTT, namely, k = arccos[maj/(4Jo(2Vb/w))]. A 
solution for k = -in: 12 can be found only for m = 0. This 
behavior is identical to that obtained from the Magnus 
expansion for w 3> Vq. From Fig. 21 we find that for high 
Vo, Floquet perturbation theory works better at higher 
frequencies where Umax = 2Jo(2Uo/w). We note that for 
very small values of w, Umax falls off as I/Vq as predicted 
by the adiabatic theory. Therefore, a crossover in the 
behavior of Umax as a function of Vq and uj is expected. 
The crossover happens between two types of behaviors 


of the maximum group velocity, i.e., Umax oc and 
Umax (X 2Jo(2Uo/u;). Although Umax never becomes zero 
(but shows a dip) at the zeros of a Bessel function for 
small UJ, we find that Umax indeed vanishes at these points 
for higher values of ui. 

The vanishing of Umax when Jo(2Uo/wo) = 0, corre¬ 
sponds to the coherent destruction of tunneling^i^ or 
dynamical freezing^. When the bare energies (diagonal 
terms) of a two-level system are sinusoidally driven with 
a driving frequency which is much larger than the tun¬ 
neling (appearing in the off-diagonal terms of the corre¬ 
sponding 2x2 Hamiltonian), the system may get frozen in 
its initial state even though the dynamics is perfectly uni¬ 
tary; this is known as the coherent destruction of tunnel¬ 
ing which occurs when the transition probability to the 
other state given by Jo(2Uo/wo) vanishes. The present 
model ideally represents such a situation when Uq 1 
and Wo 1 with Uq/wq ~ 0(1) as also observed numer¬ 
ically in Fig. 21 Whenever Umax vanishes, the quantum 
correlations do not propagate as we will discuss in Sec. 
Em this also leads to a real space localization of hard 
core bosons. 
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FIG. 4: (Color online) Variation of the maximum group veloc¬ 
ity as a function of Vq. For higher values of Vq and frequency 
Lu, the numerically obtained Umax is found to match the Bessel 
function given by 2 Jq{2Vq/uj). On the other hand, for very 
small values of the frequency Unax does not match the Bessel 
function. 


VI. SMALL POTENTIAL COMPARED TO THE 
HOPPING AMPLITUDE: Uo < 7 


We now consider the other limit, Ub 7 = 1, when Vq 
can be treated perturbatively. At the zeroth order in Uq, 
the Hamiltonian m reduces to Hk = — 2 cosA: cr^ with 
eigenfunctions 


m^)) 


COS k t \ 

0 j’ 


„i2 cos k t 


(14) 
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In the non-degenerate case when ^ ^ gi2cosfc 

T = 21: juj, it can be easily shown that the first order cor¬ 
rection in the quasi-energy vanishes since (cr®) = 0 when 
the expectation values are calculated with the eigenfunc¬ 
tions in Eq. m- This necessitates the application of 
a degenerate perturbation theory^ when the condition 
g-i2cosfc T _ jg satisfied, implying that 4cosfc = muj. 
We will distinguish between two situations, 4 cos k/uj ^ 1 
and dcos/c/w = 1; as we will show below, in the former 
case there is a correction of order Vq, while in the latter 
case a correction of order Vq emerges. 

Let us first discuss the situation in which 4 cos fc/w yf 1. 
In the same spirit as in Sec. El the quasi-states are chosen 
to be 


Ww) = 

W) = ( g-(y) g'?2cosfc t ) ■ (15) 

We note that \4>ti'r)) = and 1^^(T)) = 


gi2Tcosfc|^^ time-dependent coefficients satisfy 

the Schrodinger equation 

ic^(t) = -fVbsin(a;t) c^(t) 6“*"“ 

ic^(t) = -iVbsin(a;t) c^(t) (16) 


Within the first order perturbative approximation, we 
substitute cf(t) = c^(0) on the right hand side of the 
above equations. At t = T, we find that 


4 m = 4(0)cost), 
c;{T) = c; (0) + ^ ^ ain(4rcoa t). ( 17 ) 

Considering first the situation 4cosfc w, the Floquet 
operator up to the first order in Vq at time t = T is given 

by 


Fk{T) 


^—i2T cos k 

, sin(4r cos ' 

0 ;*^ —lo cos^ k ^ ' 


- 0 -iu sin(4r cos /c)e-*2^ 1 

— lo cos-^ k ^ ' 

^i2T cos k 


(18) 


Diagonalizing the Floquet operator in Eq. m, we get 
the Floquet quasi-energies exp(f^^r) = cos(2rcosfc) ± 

i\JsiT?{2Tcask) + [a;Fosin(4T cosk)/{uj‘^ — IGcos^ fc)]^, 

and hence = 2 cos A: -|- 0{Vq). The maximum group 
velocity becomes equal to 2 for small Vq. 

On the other hand, when 4cosfc = w, we have 
{T)) = \(j)+ (0)) and \<j)l (T)) = \(j)l (0)) to zeroth order 
in Vq. Solving the Schrodinger equations within the first 
order approximation, one finds that the time-dependent 
coefficients are given by 


4m=ct(o)+dfiA, ,-(T) = cno)-^iL(®. 

(19) 

The Floquet operator is given by 


Fk{T) 


1 Vq-kIuJ 
—VqttIuj 1 


( 20 ) 


The eigenvalues of the Floquet operator are = 1 ± 

iVQ'Kjid. The quasi-energy fj,^ = log(l ± iVott juj) jT ss 
±Fo/ 2, leading to a first-order correction to the quasi¬ 
energy unlike the previous case 4 cos k ^ uj. We also find 
that the quasi-degenerate momentum modes are given 
by fc = ± arccos(mw/4). Referring to Fig. [5] (b) for 
Vo = 0.5 and w = 27r/10, we note that the number of 
quasi-degenerate points is successfully predicted by this 
theory. A correction to the quasi-energy of the order of 
Vq appears at only m = 1. 


VII. LIGHT CONE LIKE PROPAGATION OP 
PARTICLES IN REAL SPACE WITH 
STROBOSCOPIC TIME T 


In the earlier sections, we discussed the maximum 
group velocity V^ax for a given set of parameter values 
Vq and w as presented in Fig. [TJ Here we illustrate how 
the light cone effect arising due to the existence of an 
upper bound to the group velocity manifests in the real 
space propagation of particles as shown in Fig. [Sj we see 
that there is a dynamical localization when Vmax 0. 
This is illustrated by choosing an initial state at t = 0 of 
a 200-site system in which the sites labeled 51 to 150 are 
filled (shown by the light region) and the remaining sites 
are empty (shown by the dark region); this initial state 
evolves with the total Hamiltonian, i.e., the tight-binding 
part as well as the sinusoidal driving of the staggered 
potential. For every stroboscopic instant (t = nT), we 
can find the particle density at each site by numerically 
studying the time evolution of the initial density matrix 
p(0), namely, p{nT) = F{nT)p{0)F^nT), where F{nT) 
is the real space Floquet operator at time t = nT. The 
slope of the red dotted line separating the occupied and 
unoccupied regions in Fig. is proportional to ±Vmax; 
this clearly demonstrates the light cone like propagation. 

The above scenario leads to a couple of important ob¬ 
servations. First, when the system is observed strobo- 
scopically one finds a linearly spreading boundary sepa- 
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FIG. 5: (Color online) Density of particles as a function of the stroboscopic time nT and lattice site L. The red dotted line 
signifies the localization-delocalization boundary. One can determine the maximum group velocity from the slope of this line 
[Knax = dL/{Tdn)\, which is found to agree well with Vmax obtained analytically from a Bessel function. Fig. (a) depicts a 
situation where a dynamical localization (dynamical freezing) nearly happens, for oj = drr and Vb = 15; I4nax nearly vanishes 
when Jo(2Vb/cu) = 0. (b) shows that particles move with a maximum group velocity of Vbnax = 0.814 for a; = dvr and Vb = 24. 
(c) and (d) show that no dynamical localization is observed for = 27r/3 with Vb = 12 and Vb = 14, respectively. 


rating the occupied and unoccupied regions. This empha¬ 
sizes the existence of a 14nax with which information (in 
this case, the bosons themselves) can propagate, thereby 
establishing an equivalent to the Lieb-Robinson limilii in 
a sinusoidally driven quantum system which is at a gap¬ 
less quantum critical point in the absence of the driving 
term. Second, there can be situations when 14nax van¬ 
ishes in the asymptotic limit n —>■ oo, which corresponds 
to a real space dynamical localization of the particles; 
this also implies that the system stops absorbing energy 
over a complete period even if it is being periodically 
driven. However, the particles would spread uniformly if 
the driving is stopped, which leads to the conclusion that 
this localization is indeed a result of the periodic driving. 

Recently, there has been an experimental observation 
of light cone like spreading of “two-point parity corre¬ 
lation” in an optical lattice under a sudden quenching 
of the on-site interaction strength from a deep Mott in¬ 
sulating phase to the vicinity of a superfluid-Mott in¬ 
sulator boundary^. The existence of an upper bound 
on the speed has been explained using the notion of the 
counter-propagation of quasiparticles (“holon” and “dou- 
blon”) generated due to the quench. Moreover, a dynam¬ 
ical localization-to-delocalization transition has been ob¬ 
served in a quantum kicked rotator, realized by placing 
cold atoms in a pulsed, far-detuned, standing wave, by 
measuring the number of zero velocity atoms under the 
influence of a quasiperiodic driving^. In connection to 


our work, the dynamical localization we predict can be 
experimentally observed by realizing the hard core boson 
model in an optical lattice with a sinusoidally varying 
alternating on-site potential and measuring the current 
stroboscopically starting from an initial current carrying 
ground state (obtained by applying a synthetic gauge 
potential^ to the one-dimensional optical lattice). A 
vanishing current in the large time limit {t = nT with 
n 2> 1) for a particular set of values of Vq and lu would sig¬ 
nify the existence of a dynamically localized state. Sim¬ 
ilarly, the light cone propagation of a wave packet can 
also be realized by measuring the quasiparticle correla¬ 
tion function with time; the upper bound on the velocity, 
14nax, can be determined from the first maximum of the 
correlation function as a function of time and distance. 


VIII. CONCLUDING REMARKS 

We have analyzed the behavior of a one-dimensional 
system of hard core bosons which have a nearest neigh¬ 
bor hopping amplitude 7 = 1 and are driven by a sinu¬ 
soidally varying staggered potential with magnitude Vq 
and frequency lu. We have derived the maximum group 
velocity Knax from the quasi-energies computed numeri¬ 
cally from the Floquet operator. A number of analytical 
approximation methods have been used to study Umax 
























in different regions in the parameter space. Within the 
adiabatic approximation (which is valid when a; —>■ 0), 
we find that Vmax is independent of w for small Vq and 
scales as 1/Vb for large Vq. For large frequencies we use 
the Magnus expansion of the Floquet Hamiltonian and 
find that T4iax = 2, independent of the magnitude of Vq. 
In this limit, the periodic perturbation vanishes on aver¬ 
age and only the tight-binding part of the Hamiltonian 
contributes to the group velocity. In the limit Vq 3 > 7, 
we show that the Floquet perturbation theory correctly 
predicts the vanishing of I4iax for Jq{2Vo/uj) = 0; this 
dynamical localization is particularly prominent in the 
limit of large Vq and w with Vq/uj In the other limit 
Vb 7, there is a correction to the group velocity at first 
order in Vb when the condition 47 cos k = muj is satisfied. 

Finally, starting from an initial state where the hard 
core bosons are localized in one part of the chain, we 
demonstrate that the existence of I4iax sets an upper 
bound on the speed with which particles propagate lead¬ 
ing to a light cone like spreading of the particles in real 
space. The dynamically localized (or frozen) phase is 
characterized by a vanishing I4iax- 

None of the analytical methods work in the interme¬ 


diate region when Vb and w are both of the order of 7 
(shown by the central region in the right panel of Fig.[T|). 
An analysis of the behavior of Vmax in this region may 
be an interesting subject for future research. 

We conclude with the remark that the result presented 
here for a sinusoidal driving is not special to a one¬ 
dimensional model of hard core bosons (which is equiv¬ 
alent to a fermionic model in one dimension). A similar 
behavior of the stroboscopic group velocity, especially the 
dynamical localization, can be shown to occur in mod¬ 
els with non-interacting fermions on a variety of higher 
dimensional hyper-cubic lattices. 
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